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O ■ An unambiguous identification of glueballs in experiments will be of great significance, 

because their existence is an important test of QCD. The proposal, advanced here. 
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O [ is to experimentally search for glueballs as peaks in the invariant mass of a leading 

^ I i^g-pair fragmenting from an energetic gluon jet out of high-statistics three-jet events 

■^ ■ in hadronic decays of the weak neutral Z boson. Using a physically motivated model 

p ■ of the gluon-glueball fragmentation function, we find a substantial fragmentation rate 

jL . into a leading glueball. It is very likely that a search, along the lines suggested here 

D ! by any of the four groups at the Large Electron Positron collider at CERN, will prove 

^. ' fruitful. 
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Gluons, the confined colour-octet mediators of strong interactions in quantum cliro- 
modynamics (QCD), 'sliine |l[] in tlieir own light'. They have self-interactions as a 
consequence of non-abelian gauge symmetry. At short distances << lGeV~^, their 
coupling strength decreases via renormalisation group evolution to yield an asymptot- 
ically free weak-coupling description. But at distances > (Aqcd)~^, where Aqcd 
is the QCD scale ~ 200 MeV, the coupling strength increases to a strong enough 
value to cause colour confinement. With such strong couplings, colour-singlet gluonic 
bound states or glueballs |l]] are expected to form. Indeed, there exist strong theoretical 
arguments ^ favouring such formation. 

Simple representations of scalar, pseudoscalar and tensor glueball fields are G{x) ~ 
TrFf,^{x)F^"'{x), G{x) ~ TrF^^(x)F'^^(x), and G^^{x) ~ TtF^p{x)FP''{x), respectively. 
Here F^^{x) is the covariant colour-contracted gluon field-strength tensor in standard 
notation and F^'^^x) is its dual. Theoretical studies, carried out over two decades, 
suggest the existence of glueball states in the few GeV mass-range. Such considerations 
cover bag 0] , quasiparticle and instanton models and even those on supergravity 
1^. There have also been glueball simulations [Q on the lattice which fall in the same 
ambit. 

All of the above studies predict the lightest member of the glueball spectrum to be 
a scalar. In fact, the lattice approach P] pins down its mass within the window 1.5-1.7 
GeV. A scalar mass of 1.5 GeV or so is suggested anyway from the square- root of 
the inverse of the slope ~ 0.4/GeV^ of the Pomeron trajectory, describing high-energy 
diffraction, if the latter is identified as the grandparent of the scalar glueball trajectory. 
However, there is controversy over the predicted spin, parity of the lightest glueball. 
Estimates [^ from QCD sum-rules show preference for the latter being a pseudo-scalar, 
while some field-theoretic models []TU| suggest that it could be tensor. Glueballs will. 



of course, be unstable against hadronic decay. But, on account of the VOZI rule |]l[ 
their width should not be much more than 100 MeV or so. Thus they are expected to 
be narrower than typical qq resonances in that mass range, though this characteristic 
feature may get diluted due to glueball-meson mixing. 

Since glueballs are inherently quantum chromodynamic in nature, the confirmation 
of a glueball would constitute direct evidence for QCD. Much effort, as reviewed in 
Refs. []I], |Tl|, |T2|, has gone into the production and detection of such states. First of 



all, glueballs need glue-rich production channels. They are scarcely produced in usual 
quark-antiquark creation, annihilation and rearrangement subprocesses. A further 



complication is that |13] glueballs are expected to mix significantly with flavour- singlet 
qq mesons of the same spin and parity. For instance, the central region of hadropro- 
duction is characterised by the gg production channel. Nevertheless, careful filtering 



procedures ||14[ have to be devised to avoid misidentifying flavour singlet mesons as 
glueballs among resonances produced here. There have been several quests in this di- 
rection [|l^. Another probe |l^, |l^ has been the radiative decay of charmonium {J/ip) 



where the photon could recoil against a glueball. However, the reduced statistics of 



a radiative process constitute a limiting factor here. Several interesting candidates 
have emerged from both of these studies : /o(1500), /j(1710) with J = 0,2, ^(2230) 
etc. The glueball interpretation of these flavour singlet mesons is quite plausible. Still, 
statistically significant clinching evidence for a conclusive glueball identification has 
been lacking so far and alternative avenues need to be explored. This motivates us to 
propose a new way of gathering such evidence. 




Figure 1: Glueball fragmentation from a hard gluon in Z — > qqg decay 



We are guided by the simple ides that a sufficiently hard gluon, hadronizing as 
an energetic jet, will naturally fragment first into a leading glueball. Because of the 
reasonably large mass of the glueball, a sizable amount of rapidity will be taken up 
in the process of lifting a glueball from the vacuum, leaving the residual gluon as 
quite soft (Fig. 1). Such a possibility was first mooted [|1^ two decades ago in the 
context of coUinear hadronic decays of a heavy quarkonium. Subsequent experimental 
searches on the T resonance were unsuccessful, but then the three gluons emerging 
from the bound hb annihilation in hadronic T decay are not sufficiently energetic to 
form isolated jets, which are necessary for the fragmentation process. With higher 
energy gluon jets, fragmentation processes become important as is evidenced by the 
study of J /ip production via gluon fragmentation at the Tevatron |TP[. In fact, at high 
energy colliders gluon fragmentation becomes dominant and becomes an important 
discovery channel for new particles. We are thus naturally led to direct our attention 



to hard, isolated gluon jets in the sample oi Z -^ qqg three-jet events at LEP. The 
least energetic of the three jets in the sample is taken, with a high degree of reliability 
(with an efficiency of about 70%), to be a gluon jet. Indeed, this is what is borne out 
in simulations pO[ based on perturbative QCD. Even after the imposition of a cut of 



ii^jet > 15 GeV in the Z rest-frame, one should still be left with nearly hundred thousand 
events from the LEPl Z sample. This is a rich repository of events containing isolated 
hard gluon jets. One is likely to get an observational handle on any glueball produced 
in them if one can estimate the gluon-glueball fragmentation probability, multiplied 
by the branching ratio for the glueball decaying into a i^5-pair, that is credible even 
within an order of magnitude. The glueball will show up as a peak in the KsKs 
invariant mass spectrum, studied in the gluon (but not in the q— or q—) jet. For a 
scalar or pseudoscalar glueball, no correlations are expected between the /Tg-directions 
and the jet axis; for a tensor one there will, in general, be such correlations. 

Let us quantitatively consider the question of the glueball fragmentation of a hard 
gluon in the three-jet final state of Z hadronic decay at LEP. We will consider the 
quantity V{Z -^ qqGX), the partial width for the Z to decay into a qq pair plus a 
glueball state, G, and other soft gluons. These soft gluons will produce soft hadrons 
which in this inclusive mode we denote by X. The hardness of the fragmenting gluon 
can be ensured by a cut on the energy of the gluon jet in the rest frame of the Z. It is 
expedient to write T{Z -^ qqGX), in terms of the qq partial width, T{Z — >• qq). This 
can be done using the well-known expressions for T{Z ^ qqg) in terms oiT{Z -^ qq), 
given as 

dT{Z -^ qqg) _ 2a. -^.^ . xj + xj 

dx2dx3 3tt (1 — a;i)(l — X2)' 

where Xi = 2Ei/Mz with Ei = 1, 2, 3 denoting the energy (in the Z rest frame) of the 
antiquark, quark and the gluon jet, respectively. We note that xi + X2 + X3 = 2. Using 
this expression, we can write down the corresponding expression for T{Z ^ qqGX). 
It is more convenient to write this width in terms of z, which is the fraction of the 
parent gluon energy carried by the glueball, rather than in terms of x^. After this 
transformation of variables, we fold T{Z -^ qqg) with the fragmentation function 
D{z, Q^), where Q^ is the scale at which the fragmentation function is evaluated. The 
resultant expression is 

T{Z -. qqGX) = ^-^T{Z -. qq) j dx, f dz^ . ""'^.f ^ D{z,Q\ (2) 
6Ti J J z (1 — xij^l — a;2) 

where the limits of integration in the above expression are chosen in a way consistent 
with the experimental cuts to be specified in detail below. 

To estimate the glueball production rate, we need to make an ansatz for the glueball 
fragmentation function. The simplest assumption is to consider the fragmentation of a 
high energy gluon into a glueball as being analogous to the fragmentation of a valence 
quark into a meson. This may appear unusual at first sight. We know that all quark 



jets predominantly fragment into mesons whereas in most of the gluon jets - studied in 
three-jet samples in e~^e~ machines at CM energies far below the Z-mass - the parent 
gluon first goes into a qq pair which hadronize in terms of vr's, r/'s, p's etc. Unlike 
the former, which is a zeroth order process, the latter is 0{as) in the rate; but the 
large mass of the glueball makes it impossible for such a gluon to effect a zeroth order 
fragmentation into it. Our claim is that, once a gluon is very energetic, as is the case for 
the one emitted by the Z via Z —>■ qqg, it will easily overcome this threshold effect; its 
fragmentation into a glueball state would then become a 'valence-like' process. For such 
a gluon, fragmentation via the transition first into a qq pair would be comparatively 
down by an 0{as) factor just as quark fragmentation via gluon radiation is smaller as 
compared with the direct fragmentation into a meson of a valence quark. 

There is more justification for the above assumption. A calculation |^, based 
on QCD sumrules, of the exclusive distribution of gluons inside a glueball (i.e. the 
wavefunction) shows that the results are very similar to that of the meson wavefunction. 
In fact, these calculations suggest a somewhat larger normalisation for the glueball 
wavefunction and it is possibly true of the inclusive fragmentation function too. But 
without dwelling too much on these finer points, let us point out that this ansatz for the 
fragmentation function is being made with the idea of estimating the rate of glueball 
production via fragmentation at LEP2 energies. The assumption we make allows us 
to make a rough estimate for the number of glueball events we expect to see at LEP2. 
Taking this number as given, we can then try to understand whether it is feasible to 
attempt a search for the glueball state through its decay into mesons. In that sense, 
we should take the numerical results presented here as a rough guide to decide on what 
kind of search strategies will be appropriate in the experimental situation. Moreover, 
we also present results with a different fragmentation function and study the effect 
on our results of varying this input. For the pion fragmentation function, we use 
the parametrisation of Ref. [^ (which is a 1 — z distribution with the normalisation 
obtained from a fit to pion production data) and use this as the glueball fragmentation 
function at the input scale /io ~ 2 GeV. To take into account the fact that the glueball 
mass is quite substantial, we multiply this fragmentation function with a multiplicative 
threshold factor {1 — 4Mq/E^), where Mq is the mass of the glueball and Eg is the lab- 
frame gluon energy. The fragmentation function is then evolved to the scale typical 
of the fragmenting gluon using Altarelli-Parisi evolution. In the evolution, we have 
neglected the non-diagonal anomalous dimensions, since their effects are sub-leading. 
When we vary the fragmentation function, we choose a (1 — z)^ distribution instead, 
but we normalise the distribution in such a way that the integrated probability is the 
same as in the case of the 1 — z distribution. 

In order to make contact with the experimental jet selection criteria used in the 
LEP experiments, we require that the lowest energy parton is identified as the gluon 
and that is the fragmenting parton. Also, it is usual to select the jet sample by requiring 



a ininimum cut, (imin on the quantities dij, defined as 

2EiEjSmeij/2 

Ei + Ej 



dij 



(3) 



a: 



where the i, j indices refer to the three partons in the three-jet final state. Following 
the experimental cuts, we take dmin to be 7 GeV. In addition, we also require that the 
gluon energy be above a minimum value, £'cut- Since V{Z ^> qqGX) is a function of 
-Ecut, we study this functional dependence by varying E^nt- 



0.0014 



0.0012 



0.001 



S 0.0008 



0.0006 



0.0004 



0.0002 




E cut 



Figure 2: Branching ratio into qqG final state times branching fraction of the G 
KgKg decay 



Because of the good reconstruction efficiency for the Ks at LEP, we focus on 
the decay of the glueball into KsKs, rather than for its decay into r^'s for which the 
efficiency is rather poor. Theoretical estimates |23| for the decay branching ratio of 
the glueball in the KsKs channel suggest that this could be conservatively placed at 
about 2.5%. We present our results in terms of the branching ratio into the qqG final 
state, with the branching ratio of the glueball decay into KsKs also folded in. Thus 



we define 

_ r(Z ^ qqGX) T{G ^ KsKs) 

T{Z -^ qq) V{G -. all) ' ^ ' 

In Fig. 2, we have shown our results for this branching ratio as a function of the cut 
on the energy of the gluon jet denoted as -Ecut- These are shown for both sets of 
input fragmentation functions - the curve marked Set I is with the pion fragmentation 
function and that marked Set II is the (1 — zY fragmentation function. Assuming 
four million hadronic Z's and folding in a Ks reconstruction efficiency factor (which is 
taken to be 18%), we find that with a Ec^t of about 15 GeV one would expect of the 
order of 100 events in the KsKs channel, for the Set I fragmentation function. For the 
Set II fragmentation function, this number varies by about 10%. Thus it should be 
possible for any of the four LEP groups to mount a glueball search on their three-jet 
hadronic events from the Z. As mentioned earlier, we have made a rather conservative 
choice for the normalisation of the fragmentation function. If the normalisation of the 
exclusive distribution amplitude for glueballs relative to that of the pion |2T[] is taken 



as a bench-mark, then we could expect a larger normalisation for the fragmentation 
function and a correspondingly larger number of glueball events in the hadronic decay 
of the Z. 

We also find, from our computations, that the z values that are sampled in the 
fragmentation process lie in a not very broad range at relatively small z between 0.05 
and 0.25. The lower value of z accessed, is close to the kinematic lower limit. The 
existence of the upper cut-off on z suggests that it may be able to improve the efficiency 
of the glueball search by restricting the lab energies of the glueball to be less than 
about a quarter of the energy in the gluon jet. As mentioned earlier, we had used a 
multiplicative threshold factor in the fragmentation function. But we find that the cuts 
on dij and on the gluon energy ensure that the energies involved in the fragmentation 
process are large enough, so that the effect of this factor is negligible. A similar 
kinematic behaviour results in the production of quarkonia through fragmentation of 
high energy gluons ||19| . 



We would like to emphasise that the search we have proposed in this paper is for a 
pure glueball state. It is, however, quite likely that a glueball state in the mass range of 
1.5 or 2 GeV may mix with scalar isosinglet qq states in the same mass region. Indeed, 
such a mixing has been invoked in the analysis of the /o(1500) - a glueball candidate. 



It has been pointed out |2J] that the mixing of the scalar glueball state with a qq state 
nearly degenerate in mass can change the glueball couplings so that the decays of the 
mixed state need not be such as to give equal fraction of vr's and K^s, as would be 
expected in the case of the decays of the pure glueball state Q In the event that the 
mixing is substantial, we would expect that the KgKs branching ratio of the mixed 

^The expectation that a pure glueball state would decay into equal fractions of tt's and K^s is 
rather naive and ignores the decay dynamics. The lattice calculation in Ref. [ p5| , in fact, shows that 
these rates are unequal and, in particular, in agreement with the experimental data on /j(1710). 



state to be reduced from the value used in the present calculation by the square of the 
cosine of the mixing angle. 

We have, in this letter, proposed a new way of exploring a glueball in the fragmenta- 
tion of hard gluon jets at LEP. Our estimated numbers do look sufficiently encouraging 
for any of the four LEP experiments to mount a glueball search in this channel. The 
observation of the glueball state will provide a confirmation of one of the important 
non-perturbative predictions of quantum chromo dynamics. 
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